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Abstract 

We explore the ability of the recently established quasilocal density functional the- 
ory for describing the isoscalar giant monopole resonance. Within this theory we use 
the scaling approach and perform constrained calculations for obtaining the cubic and 
inverse energy weighted moments (sum rules) of the RPA strength. The meaning of 
the sum rule approach in this case is discussed. Numerical calculations are carried 
out using Gogny forces and an excellent agreement is found with HF + RPA results 
previously reported in literature. The nuclear matter compression modulus predicted 
in our model lies in the range 210-230 MeV which agrees with earlier findings. The 
information provided by the sum rule approach in the case of nuclei near the neutron 
drip line is also discussed. 



PACS numbers: 21.60Jz, 31.15Ew, 31.15Gy 
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1 INTRODUCTION 



Recently we have established the quasilocal density functional theory (QLDFT) and its 
application for describing the nuclear ground state properties II] . It is just a generalization 
of the local Hohenberg-Kohn-Sham (HKS) theory El and is based on the definition of 
an universal energy density functional £ of the Slater density matrix (DM) po [Xj- One 
can define the quasilocal energy functional E[h, f, J] through a many-to-one mapping of the 
density matrix p to a set of local quantities n = {n p , n n }, f = {r p , r n }, J = { J p , J n }, where 
n q , T q , J q are the local, the kinetic energy, and the spin densities of each kind of nucleons 
(q=p,n): 

E[h,r,J] = inf _£[p ]. (1) 

The main property of the functional {1} is that its minimum provides the exact ground state 
energy, which is attained at the true local nucleon densities n. It is worthwhile to note that 
in this case the equilibrium f and J densities are not the exact ones and they correspond to 
the system without correlations. 

The functional E[n,f, J] consists of two terms: 

E[h,t,J) = E [n,t,J] + E RC [n], (2) 

where E Q [n, f, J] is a Hartree-Fock (HF) energy functional reduced to the quasilocal form, 
and Erq is the residual correlation energy. The HF contribution corresponds to a finite- 
range density-independent effective force. Its quasilocal form is calculated using the extended 
Thomas- Fermi theory [3], and consists of a kinetic energy part, a Hartree term, a local 
exchange contribution, and a spin-orbit energy as the used in the Skyrme energy functionals 
[T]. To compute this HF functional we use the density- independent finite-range part of the 
Gogny force The main formulas for Eq are reported in pp. The residual correlation 
energy Enc[n] is taken phenomenologically and is parametrized similar to the contribution 
of the density-dependent part of the Skyrme or the Gogny interactions. 

The aim of the QLDFT is to describe properties of the nuclear ground state such as the 
total energy E GS , the local particle densities n, the neutron and proton separation energies 
(chemical potentials) of double magic nuclei. To do that use is made of the variational 
principle to obtain a set of single-particle equations similar to the Kohn-Sham ones but 
containing additionally a position-dependent effective mass and a spin-orbit potential pQ. 
We have used this quasilocal approach to obtain some nuclear ground-state properties of 
several magic nuclei with the Gogny D1S force jH]. We have found a very good agreement 
between the binding energies and root mean square radii computed using the QLDFT and 
the corresponding values calculated with the full HF method pQ. 

Now we are in turn to investigate if the QLDFT is also able to predict some informa- 
tion about excited nuclear states. In this paper we want to analyze the collective breathing 
mode and to study the isoscalar giant monopole resonance (ISGMR) in finite nuclei using our 
model. The experimental value of the excitation energy of the ISGMR is mainly extracted 
from the analysis of the inelastic (a, a') scattering data (see [7] and references therein). From 
the experimental excitation energies in medium and heavy nuclei it is possible to estimate 
the nuclear matter compression modulus K^. In particular, self-consistent HF plus random 
phase approximation (RPA) calculations using Gogny |H] or Skyrme forces |S] determine 
to be 210-230 MeV. 
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2 THE BASIC THEORY 



Giant resonances are understood in terms of small amplitude oscillations of nuclei as a 
response to an external field generated by electromagnetic or hadronic probes. The most 
widely used theoretical framework for describing these vibrations is the RPA ^Hj which allows 
to obtain the strength function S(E) that measures the nuclear response. For medium and 
heavy nuclei far from the drip line, the strength S(E) corresponding to the breathing mode 
is mainly concentrated in a rather narrow region of the energy spectrum. In these cases the 
knowledge of a few low energy moments of S(E) (sum rules), which are defined as 

m k = J2(Eu-E ) k \(u\Q\0)\\ (3) 

where Q = ^ rf is the one-body monopole excitation operator, can provide a useful infor- 
mation on the average properties of the ISGMR. If k is an odd integer, the sum rules 
can be written as expectation values of some commutators calculated in the exact ground 
state [TT]. For example 

mi =<0\[Q,[H,Q]]\0> (4) 

and 

m 3 =<0\[[Q,H],[H,[H,Q]]]\0> ■ (5) 
From these moments one can estimate several average energies as: 



Ek,k-2 = \ — — ■ (6) 

V m k-2 

The full quantal calculation of the sum rules is still a complicated task because the 
exact ground state wave function is usually unknown. However, if the energy moments of 
S(E) are evaluated at lplh (RPA) level, it is possible to replace the exact ground-state wave 
functions by the uncorrelated HF ones in the calculation of the sum rules. 

For forces which commutes with the monopole excitation operator Q, as it happens 
for the Skyrme and Gogny interactions, the only contribution to the commutator [Q, [H, Q}} 
comes from the kinetic energy part of the Hamiltonian. Thus the mi sum rule is given by 
(seeHH) 

2h 2 „ 

m 1 = —A<r 2 > , (7) 

where M is the nucleon mass, the expectation value of the operator r 2 is calculated with the 
HF wave functions. 

The direct evaluation of the commutators entering the formula © for 7773 is a rather 
cumbersome task which can be avoided computing the 777.3 sum rule of the RPA monopole 
strength function through the scaling method [TT] . Using the scaled ground-state wave 
function 

*,(ri, ...,r A ) = J7W fori, . . . , rjr A ) , (8) 

the plus three energy moment can be expressed by means of the second derivative of the 
scaled ground-state energy Eq S {t]) =< *ff v \H\ty v > as 

/ TJ=1 



2 V M J V d V 
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Thus the 7713 moment measures the change of energy of the nucleus when the ground-state 
wave function is deformed according to (JHJ). 

With account of Eqs. (0), ©, an d Q, the scaled average energy of the ISGMR is 
defined as: 

E h MR = ^- (io) 



This method allows also to define the scaled nuclear compression modulus for a finite nucleus 
with mass number A as (see [T2] ) 

ft? >o [E GMR ) - A y ^ 2 . (11) 

Let us describe the method of calculation of the derivatives in the right-hand sides of 
Eqs. (jHl) and Within the QLDFT the scaled energy can be written as 

E S GS (V) = T( V ) + E^M + Eflgfr) + E Coul (v) + E so (v) + E RC (v)- (12) 

Due to the fact that in the QLDFT we deal with Slater determinant wave functions, the 
different contributions to the scaled energy can be easily determined. The kinetic, the spin- 
orbit, and the Coulomb contributions scale as: 

T(r ] )=r ] 2 T(l), E so {rj) = rfE so {l), and E Coul ( V ) = V E Coul (l) . (13) 

In our approach the residual correlation energy is chosen phenomenologically (see Introduc- 
tion) as: 

ER C [n] = t j J drn a (r)[(2 + x 3 )n 2 (r)-(2x 3 + l)(n 2 p (r)+n 2 n (r))}, (14) 

which under the scaling transformation reads: 

E RC {i 1 )=^ a+l) E RC {l). (15) 

These formulas give a simple way to obtain some of the terms entering ()12j) . Sim- 
ilar explicit expressions for the direct and exchange energy contributions coming from the 
finite range part of the effective force cannot be derived. However, one can calculate their 
corresponding derivatives with respect to the dimensionless collective coordinate r] using the 
following trick. For a sake of simplicity we will consider a simple Wigner central force. In 
our approach we use Gogny-like forces with Gaussian form factors v(r/fi) where \x is the 
range of the force. In this case the direct part of the energy coming from the finite range 
effective force has the general form 

E%Li = \f drdr'n(r)v(\r - r'\/^)n(r f ) . (16) 

It is easy to see that the scaling of the direct energy is reduced to renormalization of the 
range /1: 

E§UV) = \\ drdr'n(r)v(\r - r'\/( W ))n(r>) = E^M ■ (17) 
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It enables us to write (see [131 El) 



drf 



— I drdr'n(r) 



2s 



dv(s/p) i 2 d 2 v(s/p) 



ds 



+ s £ 



ds 2 



n(r') 



where s 
relation: 



r — r' is the relative coordinate. In our calculations we, however, used another 



J2 rpDir 
u ^Nucl 



(v) 



drf 1 



2 " ^Nucl \P J J 



d 2 E 



77=1 



dp 1 '- 



(19) 



which follows immediately from Eq. (|17j). 

The same is true for the exchange Fock energy E^^ijf) coming from the finite range 
effective force. Our quasilocal exchange energy is obtained by replacing the exact Slater 
single-particle density matrix by the quasi-classical DM within the extended Thomas-Fermi 
approximation (ETF) pEj. It is easy to check that the quasiclassical DM Petf^ ,r') = 
Petf(R, s ) (with R — (r + r')/2) satisfies the correct scaling transformation law: 



Petf(R, s, 7]) = r] 3 p ETF (i]R, r]s) . 
Thus, because in the QLDFT the exchange energy has the general form 



we have 



and 



E 



Exch 
Nucl 



E 



(V) 



Exch 
Nucl 



~ J dRdsp 2 ETF (R, s)v(s/p) 



dRdsp 2 ETF (R, s)v(s/ (rip)) 



TpExch 
^Nucl 



(m) 



d 2 E§ x u c c 1( V ) 



drf 



d 2 E^(p') 



dp 1 '- 



(20) 



(21) 



(22) 



(23) 



The derivatives in Eqs. (|T9"|) . (|2^)l were calculated numerically. 

In light nuclei the strength function is much more spread and fragmented than in 
heavy nuclei (see for instance J3j). Therefore, to get more insight about the spreading of 
S(E) we will consider another estimate of the ISGMR mean energy provided by a constrained 
QLDFT calculation. It is well known (see ^1]) that the m_i RPA sum rule is half the 
ground-state polarizability with respect to the excitation operator Q of the ISGMR, i. e. 



m_i 



E 



E„ — En 



dR 2 (X) 
dX 



A=0 



d 2 E c GS W 
d\ 2 



(24) 



A=0 



where -R 2 (A) and Eq S (X) are the expectation values of the excitation (Q) and the Hamiltonian 
(H) operators evaluated with the ground-state wave function of the constrained Hamiltonian 
H c = H + XQ. Using Eqs. ©, 0, and (J2H), the constrained estimate of the average energy 
of the ISGMR is defined as 

(25) 



rpc 
GM R 



mi 
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The constrained nuclear compression modulus can also be denned in a similar way to that 
of Eq. (HU): 

M „ n9 if 9 d 2 E c r JR)\ , , 

K' A = T2 <r^> (E< OMR Y = - A U 2 —g-^ ) , (26) 

V / R=R 

where the function Eq S (R) is defined by the relation Eq S (R(X)) = Eq S (X), R = R(\ = 0). 

As it has been pointed out before, in Ref. [TT] was proved that the exact ground-state 
wave function can be replaced by the HF one in the self-consistent HF + RPA (lplh) sum 
rule calculation. Analogous considerations allow to state that the m 3 , m 1; and m_x sum 
rules calculated within QLDFT (i.e. using Slater determinants built up of the single-particle 
wave functions obtained from the minimization of the QLDFT energy functional), coincide 
with the QLDFT-based self-consistent RPA results. 

Finally, let us note that the RPA moments fulfill a/ m 3 /mi > mi/m > ^mi/m_i . 
Therefore, the scaled (Eq MR ) and the constrained {Eq MR ) estimates of the average energy 
of the resonance give an upper and a lower bound of the mean energy of the ISGMR mi/m®. 
The total RPA width a of the strength distribution can also be estimated from the E S GMR 
and Eq MR energies: 

1 



^<^\/(E s G mr) 2 -(E c G mr) 2 - (27) 



3 NUMERICAL RESULTS 



We have calculated the scaled and the constrained estimates of the average energies of the 
ISGMR as well as the corresponding compression moduli in finite nuclei using the QLDFT 
with the Dl' [H] and the D1S [Oj Gogny forces. In order to test our method we compare in 
Table ^ our QLDFT results with the calculations performed in Ref. |S] in the framework of 
the HF and HF + RPA approaches with the same forces. It can be seen that the agreement 

Table 1: The ISGMR average energies (in MeV) in 208 Pb computed in the framework of the 
scaled (Eq MR ) and the constrained (Eq MR ) approaches with the Dl' and the D1S Gogny 
forces using the QLDFT compared with the average energies obtained in Ref. [H] using the 
HF and the HF + RPA methods. 







QLDFT 


HF 


HF + RPA 




TpS 

GM R 


PC 

GMR 


GMR 


Ei t -i E 3 i 


Dl' 


14.51 


13.93 


14.05 


14.15 15.33 


D1S 


13.65 


13.05 


13.22 


13.34 14.16 



between the constrained QLDFT and HF energies Eq MR , and the average HF + RPA en- 
ergies Ei _i is fairly well. The agreement between the scaled QLDFT energies Eq MR and 
the values of E^\ obtained in HF + RPA approach is reasonable, though it is worse than 
for the constrained energies. From these comparisons we conclude that using our sum rule 
approach based on the QLDFT, the corresponding average energies can be confidently used 
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to theoretically estimate the ISGMR energies with Gogny forces, at least for nuclei for which 
the monopole strength shows a well defined narrow peak. 

In Table 121 we display results obtained for the nuclei 40 Ca, 90 Zr, and 208 Pb, for which 
the experimental excitation energies of the ISGMR are accurately known F° r the 

Table 2: The compression moduli K A (in MeV) and the ISGMR average energies Eqmr 
(in MeV) of some magic nuclei computed in the framework of the scaled (K A , Eq MR ) and 
the constrained (K A , Eq MR ) approaches with the Dl' and the D1S Gogny forces using the 
QLDFT compared with the same quantities computed with the SkM* and the SIII Skyrme 
interactions using the HF approximation. Experimental average energies are taken from 
Refs. [3EEH]. 







16Q 


28Q 


40 Ca 


90 Zr 


208p b 


K'a 


Dl' 


132 


106 


146 


154 


152 




D1S 


126 


100 


137 


141 


137 




SkM* 


124 


103 


138 


146 


143 




SIII 


199 


166 


225 


243 


243 


GMR 


Dl' 


28.1 


20.1 


23.2 


19.0 


14.5 




D1S 


27.1 


19.3 


22.2 


18.0 


13.7 




SkM* 


26.8 


19.6 


22.2 


18.3 


13.9 




SIII 


34.5 


25.3 


28.5 


23.4 


17.9 




Dl' 


102 


14 


127 


144 


140 




D1S 


100 


14 


121 


132 


125 




SkM* 


94 


28 


119 


136 


133 




SIII 


132 


17 


178 


216 


214 


PC 

^GMR 


Dl' 


24.7 


7.2 


21.6 


18.4 


13.9 




D1S 


24.1 


7.3 


20.9 


17.5 


13.1 




SkM* 


23.2 


10.2 


20.6 


17.6 


13.3 




SIII 


28.1 


8.1 


25.4 


22.1 


16.8 



E e n Zn 19.2 ± 0.4 17.9 ± 0.2 14.2 ± 0.3 



nucleus 208 Pb, the QLDFT scaled and constrained predictions of the ISGMR excitation 
energy computed with the Dl' Gogny force agree very well with the experimental centroid 
value (14.2 MeV). On the other hand the QLDFT results obtained with the D1S force 
underestimate the experimental value by a 3.5 % (Eq MR ) and 8 % (Eq MR ). The reason for 
this discrepancy lies on the fact that the nuclear matter compression modulus for the Dl' 
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force (Koo = 228 MeV) is larger than that for the D1S force (K^ = 209 MeV). The situation 
is inverse for the 90 Zr nucleus where the theoretical QLDFT estimates of the ISGMR energy 
for D1S force are in a better agreement with the experimental centroid value (17.9 MeV) than 
that for the Dl' force. Thus extracted from the sum rule approach using the results for 
the nucleus 208 Pb is around 230 MeV, while the results for the nucleus 90 Zr lead to the value 
around 210 MeV, which corresponds to the the upper and lower bounds of the compression 
modulus predicted in Ref. jH]. For the 40 Ca nucleus, the ISGMR energies obtained with 
our QLDFT plus sum rule approach clearly overestimate the experimental centroid value 
both for the Dl' and for the D1S forces. In particular, for the D1S force the calculated 
scaled energy overestimates the (m^/mi) 1 ^ 2 experimental value ^Hj by a 7 % whereas the 
constrained energy is larger than the experimental (mi/m^i) 1 / 2 by a 21 %. For this nucleus 
our QLDFT model predicts a total RPA width (Eq. of 3.7 MeV. This means that 

the monopole strength is quite spread or fragmented and thus the sum rule estimate of the 
ISGMR for 40 Ca should be considered only as qualitative. 

The use of radioactive beams allows to obtain nuclei beyond the limits of /^-stability 
and to reach the neutron drip line for nuclei with Z < 8. Nuclei near the drip line are 
characterized by the small energy of the last bound levels and by their large asymmetry. It 
is expected that the properties of such nuclei, in particular the ones concerning collective 
excitations, may considerably differ from the corresponding properties of stable nuclei. The- 
oretical HF + RPA calculations carried out by Hamamoto, Sagawa, and Zhang [THj using 
the SkM* force, have shown that in nuclei near the drip line the distribution of the monopole 
strength is much affected by the presence of the low-energy threshold due to the tiny bound 
nucleons. This pattern clearly differs from the one shown by stable nuclei where the strength 
is concentrated in a single peak at least for medium and heavy nuclei. 

In particular the aforementioned authors have analyzed the monopole strength of the 
28 nucleus ^7]. We have also considered this nucleus in our QLDFT plus sum rule approach 
and their scaled and constrained energies are also reported in the Table El Although the 
sum rule approach is unable to deal with fine details of the monopole strength, it provides 
some signatures about the global behaviour of the monopole strength in nuclei near the drip 
line. In the nucleus 28 0, the scaled and the constrained energies are clearly reduced with 
respect to the corresponding values in the stable 16 O nucleus, which are also given in the 
Table 121 This reduction is particularly significant for the constrained energy and produces a 
noticeable enhancement of the width of the ISGMR. These facts point out the importance 
of the low-energy part of the monopole strength in drip line nuclei. 

In order to do a complementary check of our QLDFT estimates of the average ex- 
citation energies of the ISGMR using the Dl' and the D1S forces, we have repeated the 
sum rule analysis using the Skyrme force SkM*, which has 7^00=217 MeV slightly larger 
than the corresponding D1S value (209 MeV). The scaled energies with Skyrme forces are 
easily calculated [TT] and the constrained calculations are performed in the way discussed 
previously. The excitation energies as well as the finite nuclei incompressibilities obtained 
with SkM* are also collected in the Tabled It can be seen that D1S and SkM* practically 
predict the same values for the scaled and constrained energies. Actually the excitation 
energies obtained with SkM* are slightly larger than the ones computed with D1S, at least 
for medium and heavy nuclei where the sum rule approach can be regarded more confidently, 
due to the fact that is a little bit larger in the former than in the latter force. In the 
same Tabled we also report the ISGMR excitation energies obtained using the Skyrme SIII 



S 



force which has 7^00=355 MeV. As it has been pointed out in [S], the energy of the breathing 
mode is roughly proportional to \fK~^,. We have checked this behaviour with the excitation 
energies of the ISGMR of the nuclei 90 Zr and 208 Pb calculated with the Dl', D1S, SkM*, and 
SIII forces using the scaling method and performing constrained calculations. It is found 
that this proportionality is also reasonably well fulfilled for each nucleus when the QLDFT 
values are used, especially in the case of the scaled energies. 

4 SUMMARY 

We have analyzed the ability of the recently established quasilocal density functional theory 
for describing the breathing mode in finite nuclei. First of all, it has to be pointed out that 
the lplh RPA calculations based on this theory have similar properties to the HF + RPA 
calculations. In particular, the exact RPA ground-state wave function can be replaced by 
the Slater determinant built up of the single-particle QLDFT wave functions in order to 
obtain the odd moments of the lplh RPA strength function. Thus the scaling approach and 
constrained calculations can be used to obtain the m 3 and m_i moments (sum rules) of the 
RPA strength based on the quasilocal density functional theory. 

We have applied this sum rule approach using the Gogny forces. For the nucleus 
208 Pb, we have found an excellent agreement with the HF and the HF + RPA results for the 
ISGMR reported in the literature. Using the Dl' and the D1S Gogny forces and comparing 
our theoretical estimates of the average energies of the ISGMR with the experimental values 
of some selected medium and heavy nuclei, we have found that the predicted nuclear matter 
compression modulus lies in the range 210-230 MeV which also agrees with earlier findings. 
It is also found that our estimates of the average energies of the ISGMR, in particular the 
ones obtained with the scaling method, roughly scale with the square root of the compression 
modulus in nuclear matter. 

Although the sum rule approach works basically for medium and heavy stable nuclei 
where the RPA strength is mainly concentrated in a narrow peak, we have analyzed its 
predictions in nuclei near the neutron drip line. In this case the RPA strength is broaded 
and the low energy part considerably enhanced because of the last weakly bound neutrons. 
The sum rule approach is able to give globally this tendency in such nuclei by reducing the 
m 3 and especially the m_i sum rules and thereby increasing the estimate of the resonance 
width. 
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